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Abstract 

Development of yellow mustard {Sinapis alba L.) with superior quality traits (low erucic and linolenic acid contents, and low 
glucosinolate content) can make this species as a potential oilseed crop. We have recently isolated three inbred lines Y1 127, 
Y514 and Y1035 with low (3.8%), medium (12.3%) and high (20.8%) linolenic acid (C18:3) content, respectively, in this 
species. Inheritance studies detected two fatty acid desaturase 3 (FAD3) gene loci controlling the variation of CI 8:3 content. 
QTL mapping revealed that the two FAD3 gene loci responsible for 73.0% and 23.4% of the total variation and were located 
on the linkage groups Sal02 and Sal 1 0, respectively. The FAD3 gene on Sal02 was referred to as SalFAD3.LAl and that on 
Sail 0 as SalFAD3.LA2. The dominant and recessive alleles were designated as LA 1 and la 1 for SalFAD3.LAl, and LA 2 and la 2 for 
SalFAD3.LA2. Cloning and alignment of the coding and genomic DNA sequences revealed that the SalFAD3.LAl and 
SalFAD3.LA2 genes each contained 8 exons and 7 introns. LA' had a coding DNA sequence (CDS) of 1143 bp encoding a 
polypeptide of 380 amino acids, whereas la 1 was a loss-of-function allele due to an insertion of 584 bp in exon 3. Both LA 
and la 2 had a CDS of 1 152 bp encoding a polypeptide of 383 amino acids. Allele-specific markers for LA 1 , la', LA 2 and la 2 co- 
segregated with the CI 8:3 content in the F 2 populations and will be useful for improving fatty acid composition through 
marker assisted selection in yellow mustard breeding. 
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Introduction 

Yellow mustard (Sinapis alba L., 2n = 24) is cultivated as an 
important condiment crop. It has many desirable agronomic traits 
such as resistance to cabbage aphids [1], flea beetles [2,3] and 
blackleg diseases [4]. In addition, it is drought tolerant and 
resistant to pod shattering. Yellow mustard germplasm with canola 
quality (low erucic acid and low glucosinolate contents) was 
developed at Agriculture and Agri-Food Canada-Saskatoon 
Research Centre (AAFC-SRC) [5], which makes yellow mustard 
have the potential to become an alternative oilseed crop to canola 
B. napus, especially in semi-arid areas. 

The oil quality of canola B. napus is determined by the 
proportion of the three major unsaturated fatty acids: oleic acid 
(C18:l), linoleic acid (C18:2) and linolenic acid (C18:3). Tradi- 
tional B. napus cultivars contain 9% C18:3 of the total fatty acids 
[6] . The high level of linolenic acid in canola oil is undesirable 
since it shortens the shelf life and causes off-type flavour of the oil 
due to the three easily oxidized double bonds. A low linolenic acid 
mutant, containing 3-5% C18:3, was produced by ethyl meth- 
anesulfonate (EMS) treatment of a high C18:3 B. napus cv. Oro 



seed [7]. Current low C18:3 canola cultivars have been developed 
using this low linolenic gene source. 

Linolenic acid content is determined mainly by the embryonic 
genotype with some influence from temperature, maternal 
genotype and cytoplasm in B. napus [8-10]. QTL mapping 
identified two major QTLs, accounting for 25.2-28.8% and 52.4 — 
62.7% of the C18:3 variation, located on the linkage groups A4 
and C4, respectively, in B. napus [11,12]. It was reported that the 
low C18:3 variant resulted from mutations of FAD3 genes in B. 
napus [11-13]. The FAD3 gene on A4 harboured a C to T 
substitution in exon 7, which when translated causes the wild type 
amino acid arginine to be replaced by cysteine. The FADS gene on 
C4 contained a G to A substitution in the 5' splice site of intron 6 
in the low C18:3 B. napus line. FAD3 allele-specific markers based 
on the sequence variation were developed and proved to be useful 
for identification of different CI 8:3 genotypes in canola B. napus 
[11,12]. Yellow mustard accessions contain 6.9-12.4% linolenic 
acid of total fatty acids in the seed [14,15]. Recendy, inbred lines 
with high (18.5%), medium (13.8%) and low (3.8%) linolenic acid 
content, respectively, have been obtained through inbreeding of 
heterozygous open-pollinated plants in yellow mustard [16]. 
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The low linolenic acid variant (3.8%) is a valuable gene source 
for breeding canola-quality yellow mustard with high stability oil 
(high oleic and low linolenic acids) as that of canola B. napus. The 
knowledge about genetic and molecular bases of the variation in 
CI 8:3 content and development of FADS allele-specific markers 
will greatly facilitate the development of low linolenic canola- 
quality yellow mustard. The objectives of this study were: 1) to 
determine the inheritance and perform QTL mapping of the 
CI 8:3 content; and 2) to clone the FADS genes and further develop 
allele-specific markers for marker assisted selection. 

Materials and Methods 

Plant Materials 

Linolenic acid contents of the three parental lines Yl 127, Y514 
and Y1035 are shown in Table 1. Yl 127 is an S 4 inbred line 
produced by selfing of the low linolenic S 2 line Y158 for two 
generations and has a low C18:3 content (average: 3.8%). Y514 is 
the doubled haploid line SaMD3 [17] and has a medium CI 8:3 
content (average: 12.3%). Y1035 is an S 4 inbred line and has a 
high CI 8:3 content (average: 20.8%). 

The Fj seeds of the three crosses Yl 127 (low)xY1035 (high), 
Y1127 (low)xY514 (medium) and Y514 (medium) xYl 035 (high) 
were produced. To produce the B^ seeds, the Ft plants of the 
three crosses were crossed as the female with the parental line with 
a lower C18:3 content. All plants were raised under the same 
conditions in the greenhouse at AAFC-SRC. 

Regional Linkage Mapping 

Regional linkage mapping of the linolenic acid content was 
performed using intron length polymorphism (ILP) markers and 
bulked segregant analysis (BSA) [18]. A total of 1478 ILP primer 
pairs: 380 from Arabidop.sis thaliana [19] and 1098 from B. napus [20] 
were used to screen the three parental lines for polymorphic 
markers. The high bulk was made by mixing equal amount of 
DNA from 10 F 2 plants with the highest C18:3 content, while the 
low bulk was formed from 10 F 2 plants with the lowest CI 8:3 
content for each of the three crosses. The primers detecting 
polymorphic markers between the two bulks were subsequently 
used to genotype individual plants of the three F 2 populations. 
Genomic DNA was extracted from young leaves of the parental 
lines Yl 127, Y514 and Y1035, F! and F 2 plants using a modified 
sodium dodecyl sulfate method [21]. Each PCR (20 u.1) contained 
1 x standard PCR buffer (NEB), 1 U of Taq polymerase (NEB), 



0.25 |J.M forward primer, 0.25 u.M reverse primer, 100 uM each 
dNTP and 50 ng of genomic DNA in a total volume 20 |xL. The 
PCR amplification consisted of an initial denaturation at 94°C for 
5 min, 35 cycles consisting of 94°C (45 sec), 55°C (45 sec), 72°C 
(1 min) terminating with 72°C for 7 min. All PCR products were 
analyzed by electrophoresis in 2 % agarose gels in 1 x Tris-acetate- 
ethylenediaminetetraacetic acid buffer. Gels were visualized by 
staining in ethidium bromide and photographed on a digital gel 
documentation system. 

The regional linkage map of CI 8: 3 content was constructed 
using JoinMap 4.0 [22] with a minimum LOD threshold of 4.0. 
QTL analysis of CI 8:3 content was performed using the interval 
mapping method of MapQTL 6.0 [23]. A Chi-square test was 
used for evaluating the genetic model of C18:3 content in the BC! 
and F 2 populations, and the ILP markers in the F 2 populations. 

Cloning of the Coding Region of the FAD3 Gene 

Primer pair No 1 (Table SI) was designed based on the 
conserved coding regions of the FADS genes in B. napus and A. 
thaliana. It was used to clone the coding DNA sequence (CDS) of 
the FAD3 gene in yellow mustard. Immature seeds at 22 days after 
pollination were collected from two individual plants from each of 
the parental lines. Total RNA was extracted from the immature 
seeds using the RNeasy Plant Mini Kit (Qiagen) as per the 
manufacturer's instructions. 750 ng of RNA from each of the 
parental lines was used to prepare the cDNA using Qiagen's 
Omniscript RT Kit as per the manufacturer's instructions. Each 
PCR (20 ul) contained 1 x PCR standard buffer (NEB), 100 uM of 
each dNTP, 0.25 |J,M of each forward and reverse primer, 1 U of 
Taq polymerase (NEB) and 50 ng of cDNA. Polymerase chain 
reaction was performed with an initial denaturation at 94°C for 
3 min followed by 35 cycles of 45 s at 94°C, 30 s at 55°C and 
1 min at 72°C with a final extension cycle of 72°C for 10 min. 

Cloning of the 5' and 3' Flanking Sequences and the 
Genomic DNA Sequences of the FAD3 Genes 

Primer pairs No 2 and 3 (Table SI) were designed based on the 
5' coding sequences of the cloned SalFAD3.LAl and SalFAD3.LA2 
genes, respectively. They were used to clone the 5' upstream 
sequences by PCR walking according to the protocol of Siebert et 
al. [24]. Primer No 4 (Table SI) was designed based on the 3' 
coding sequences of the cloned SalFAD3.LAl and SalFAD3.LA2 
genes, and was used to clone the 3' flanking sequence by PCR 
walking. Primer pairs No 5 and 6 (Table SI) were designed based 



Table 1. Linolenic acid contents of the parental lines Yl 1 27, Y514, Y1035 and seeds, and the mid-parental value in yellow 
mustard. 





Genotype 


Generation 


Linolenic Acid Content* (% of total fatty acids) 


Yl 127 


s 4 


3.8±0.7 


Y514 


DH 


12.3±0.7 


Y1035 


s 4 


20.8±0.8 


Y1127xY1035 


F, 


13.7±1.3 


Mid parent value 




12.3 


Y1127xY514 


F, 


8.9±0.7 


Mid parent value 




8.0 


Y514XY1035 


F, 


15.3±0.7 


Mid parent value 




16.5 



*: Linolenic acid content is expressed as mean value ± standard deviation. 
doi:1 0.1 371 /journal.pone.0097430.t001 
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Figure 1 . Frequency distributions of linolenic acid contents in individual seeds of Y1 1 27, Y1 035, F„ (Y11 27 xY1 035) xYl 035 and F 2 
populations. 

doi:1 0.1 371 /journal.pone.0097430.g001 
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Figure 2. Frequency distributions of linolenic acid contents in 
individual seeds of Y1 1 27, Y51 4, F n , (Y1 1 27 xY51 4) xYl 1 27 and 
F 2 populations. 

doi:1 0.1 371 /journal.pone.0097430.g002 



on the 5' flanking sequence and the 3' flanking sequence of the 
cloned SalFAD3.LAl and SalFAD3.LA2 genes, respectively, and 
were used to clone the genomic DNA sequences of SalFAD3.LAl 
and SalFAD3.LA2 genes. The standard protocol from the Clontech 
kit (website: www.clontech.com, Protocol PT 3042, Version PR 
03300) by Gwyneth Ingram and Karine Coenen was followed to 
facilitate the PCR walking. 

DNA Sequencing 

The expected PCR bands were cloned using the pGEM-T 
Vector System I (Promega) following the provided instructions. 
The plasmids were extracted using the QiaSpin Kit (Qiagen) 
following the manufacturer's instructions and sequenced using the 
primer pairs No 7-11 (Table SI) at the Plant Biotechnology 
Institute, National Research Council, Canada. 

Phylogenetic Tree 

The multiple alignments were performed using ClustalW 
(http://www.ebi.ac.uk/clustalw/). MEGA software (version 4.0) 
(http://www.megasoftware.net/index.html)[25] was used to con- 
struct a phylogenetic tree with the aligned protein sequences. The 
neighbor-joining method was used with the pairwise deletion 
option, poisson correction model, and the 1000 bootstrap 
replicates test. 
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Figure 3. Frequency distributions of linolenic acid contents in 
individual seeds of Y514, Y11035, F n , (Y514 > Y1035) Y514 and 
F 2 populations. 

doi:1 0.1 371 /journal.pone.0097430.g003 

Development of the SalFAD3.LA1 and SalFAD3.LA2 Allele- 
specific Markers 

The SalFAD3,LAl and SalFAD3.LA2 allele-specific markers were 
generated using primer pair No 12 (Table SI) which was designed 
based on the conserved flanking sequences of intron 3. The PCR 
reaction was performed with LongAmp Taq 2 X Master Mix 
(NEB) following the manufacturer's instructions with a 60°C 
annealing temperature. 

Fatty Acid Analysis 

Seed fatty acid composition was analyzed according to [26] with 
the following modification: the gas chromatography of the methyl 
esters was performed with a HP-INNOWax fused silica capillary 
column (0.25 mm by 0.5 m and 7.5 u.m) (Agilent Technologies) at 
250°C using hydrogen as the carrier gas. A minimum of 10 seeds 
from each of the parental lines and F| hybrids as well as 160 F 2 
seeds of each of the three crosses were half-seed analyzed 
according to [27]. Ninety-six seeds from each of the BCi 
populations were analyzed using the single seed method. 

Results 

Linolenic Acid Content is Controlled by Two Gene Loci in 
Yellow Mustard 

The C18:3 content of the Fj seeds was significantly higher than 
the mid-parent value in the crosses of Yl 127 (low)xY1035 (high) 
(t = 3.84,/»<0.01) and Y1127 (low)xY514 (medium) (t = 5.62, j&< 
0.01) (Table 1, Figure 1 and 2), suggesting a partial dominance of 
the high/medium over low C18:3 content. However, in the cross 
of Y514 (medium) xY1035 (high) the F, seeds had significantly 
lower C18:3 content (15.3%) than the mid-parent value of 16.5% 
(t = 6.98, /i<0.01) (Table 1, Figure 3), indicating a partial 
dominance of the medium over high C18:3 content. 



The Bd seeds of (Yl 127 xY1035) xYl 127 were classified into 
two groups: seeds with medium to high (5.4—19.1%) C18:3 content 
and seeds with low (3.0-4.9%) C18:3 content (Figure 1), fitting 
with a segregation ratio of 3:1 (% — 2.00, p - 0.16). The F 2 seeds of 
Y1127xY1035 ranged from 2.9% to 20.4% in C18:3 content 
(Figure 1) with a segregation ratio of 15:1 (seeds with 4.5-20.4% 
versus seeds with 2.9-4.3% C18:3 content) (/ 2 = 3.07, p-0.08). 
Therefore, the segregation patterns of CI 8:3 content in the B^ 
and F 2 populations supported a digenic inheritance model in this 
cross. 

The BC, seeds of (Yl 127xY514)xYl 127 showed a segregation 
ratio of 1:1 (seeds with 2.7-5.2% versus seeds with 6.4-9.7% 
CI 8:3 content) (Figure 2) (/ 2 = 3.38, p = 0.07), suggesting that the 
CI 8:3 content was controlled by one gene locus in this cross. The 
F 2 seeds of Y1127xY514 showed a continuous distribution 
ranging from 3.0% to 16.5% in the C18:3 content (Figure 2). 
The Bd seeds of (Y514xY1035)xY514 and the F 2 seeds of 
Y514xY1035 exhibited a continuous frequency distribution in the 
CI 8:3 content (Figure 3). Therefore, it was not possible to classify 
the seeds into discrete groups. 

Two QTLs Accounting for the Variation of C18:3 Content 
are Mapped to Linkage Groups Sal02 and Sail 0, 
Respectively 

In the F 2 population of Y1127 (low)xY1035 (high), eighteen 
ILP primer pairs were polymorphic between the high (16.6- 
20.4%) and low (2.9-4.0%) CI 8:3 bulks and generated 18 markers 
(Table 2). The 18 markers were mapped to two linkage groups, 
each of which carried one QTL for the C18:3 content (Figure 4). 
Based on the common ILP markers, the two linkage groups were 
revealed to be Sal02 and SallO of the constructed S. alba map [28]. 
One QTL (LOD = 45.43) accounting for 73.0% of the total 
variation of C 18:3 content was localized between BnapPIP685 and 
BnapPIP881 in Sal02 (Figure 4). The other QTL (LOD = 9.28) 
responsible for 23.4% of the total variation was located between 
BnapPIP1012 and BnapPIP363 in SallO (Figure 4). Together, the 
two QTLs explained 96.4% of the total variation for CI 8:3 
content in the F 2 population. 

In the F 2 population of Y1127 (low)xY514 (medium), 10 
polymorphic ILP primer pairs between the low (3.0-4.0%) and 
medium (14.5-16.5%) C 18:3 bulks produced 10 markers (Table 2). 
The 10 markers were all mapped to one linkage group 
corresponding to Sal02. The QTL (LOD = 46.53) was localized 
between BnapPIP685 and BnapPIP881 in the linkage group 
(Figure 4). In the F 2 population of Y514 (medium) XY1035 (high), 
1 1 markers were generated by 1 0 polymorphic primer pairs 
between the medium (10.4-11.6%) and high (16.7-19.2%) C18:3 
bulks. The 11 markers were mapped to the linkage group SallO. 
The QTL (LOD = 6.09) was located between BnapPIP1012 and 
At3g43520 in SallO (Figure 4). The two FAD3 gene loci 
controlling the QTLs in Sal02 and SallO were referred to as 
SalFAD'3 .LAI and SalFAD3.LA2, respectively. The dominant and 
recessive alleles of the SalFAD3.LAl gene were accordingly 
designated as LA and la ', while that of the SalFAD3 '.LA2 gene 
as LA and la 2 . Therefore, it could be inferred that the CI 8:3 
genotypes of Yl 1 27 (low), Y514 (medium) and Y1035 (high) were 
la'ia'la 2 la 2 , LA 1 LA 1 la la 2 and LA 1 LA 1 LA 2 LA 2 , respectively. 

The SalFAD3.LA1 and SalFAD3.LA2 Genes are Cloned and 
Exhibit Differences in the Exon and Intron 

The coding regions of the dominant alleles LA and LA were 
cloned from Y1035, while those of the recessive alleles, la 1 and la 2 
from Yl 127 using primer pair No 1 (Table SI). LA had a coding 
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Figure 4. Mapping QTLs controlling CI 8:3 content. A. The QTL in Sal02 was located between BnapPIP685 and BnapPIP881 in Y1127xY1035 
and Y1127xY514. B. The QTL in SallO was located between BnapPIP1012 and BnapPIP363 in Y1127xY1035, and between BnapPIP1012 and 
At3g43520 in Y514xY1035. 1-LOD and 2-LOD supporting intervals of each C18:3 QTL were marked by thick and thin bars, respectively. The 
SalFAD3.LA 1 and SalFAD3.LA2 genes co-localized with their C18:3 QTL peaks in the linkage groups Sal02 and SallO. 
doi:1 0.1 371 /journal.pone.0097430.g004 



DNA sequence (CDS) of 1 143 bp encoding a polypeptide of 380 
amino acids, la 1 had a CDS of 1 1 7 1 bp. Sequence alignment with 
LA 1 indicated that la harboured an indel involving a 64 bp 
insertion and a 36 bp deletion at position 412 (Figure 5 and Figure 
SI). A stop codon at the beginning of the 64 bp insertion might 
have resulted in the termination of protein translation after the 
137 th amino acid residue. Therefore, la 1 is a loss-of-function allele. 
The 5' flanking sequences from the translation start site were 



cloned for LA and la using the primer pair No 2 (Table SI). The 
5' fragment of LA was 1250 bp, while that of la 1 was 621 bp. A 
435 bp 3' flanking sequence from the translation stop codon was 
cloned for L4 ; and la 1 using the primer pair No 4 (Table SI). The 
two alleles didn't exhibit any differences in the cloned 3' flanking 
sequences. The genomic DNA sequences of the LA and la were 
amplified using the primer pair No 5 (Table SI) which was 
designed based on the 5' flanking sequence and the conserved 3' 




Exon 1 Exon 2 Exon 3 Exon 4 Exon 5 Exon 6 Exon 7 Exon 8 



Figure 5. Structure of the SalFAD3.LA 1 and SalFAD3.LA2 alleles LA 1 , la', LA 2 and la 2 in yellow mustard. The black boxes represented the 
exons and the lines between the black boxes indicated the introns. The numbers on the top of black box of LA 1 (LA 2 ) indicated the beginning and 
ending of each exon of LA' and la' (LA 2 and la 2 ). The number above each black line indicated the intron length. The nucleotide sequences of 9 bp 
deletion at position 45 of exon 1 of the alleles LA' and la', and the 64 bp insertion and the 36 bp deletion at position 412 of exon 3 of la' were 
displayed on the top. Allele specific markers were developed based on the variation in intron 3 of the SalFAD3 alleles LA', la', LA 2 and la 2 . 
doi:1 0.1 371 /journal.pone.0097430.g005 
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Figure 6. Phylogentic relationship based on the polypeptide sequences among the LA' allele of the SalFAD3.LA 1 gene, the LA 2 allele 

of the SalFAD3.LA2 gene, and the FAD3 genes of Brassica species and Arabidopsis. 
doi:1 0.1 371 /journal.pone.0097430.g006 



flanking sequence specific to the candidate SalFAD3.LAl gene. 
Comparison of the coding and genomic DNA sequences indicated 
that the candidate SalFAD3.LAl gene contained 8 exons and 7 
introns (Figure 5). Alignment of the genomic DNA sequences of 
LA and la revealed that la had an insertion of 584 bp in the 
third exon. This insertion contained a new intron splicing site GT 
(Figure S4), which resulted in a 64 bp insertion and a 36 bp 
deletion (nucleotide 412-447) at position 412 in the CDS 
(Figure 5). The inserted fragment contained a 5 bp direct repeat 
(5'-AGAAC-3') at each end, which is a typical LTR retroelement 
insertion site (Figure S4). In addition to differences in the CDS, 
LA 1 and la' exhibited variation in the length of the introns 
(Figure 5). 

Both LA 2 and la 2 had a CDS of 1 152 bp encoding a polypeptide 
of 383 amino acids (Figure SI and S5). Six point mutations at 
positions 567, 579, 666, 699, 777 and 1059 were observed in the 
CDS of la when compared with that of LA 2 , but did not lead to 
any amino acid changes. The 5' flanking sequences from the 
translation start site were cloned for LA 2 and la 2 using the primer 
pair No 3 (Table SI). The 5' flanking fragments of the two alleles 
were 444 nucleotides in length and were similar in sequence. A 
435 bp 3' flanking sequence from the translation stop codon was 
cloned for LA 2 and la 2 using the primer pair No 4 (Table SI). The 
two alleles didn't show any differences in the cloned 3' flanking 
sequences. The genomic DNA sequences of LA 2 and la 2 were 
cloned using primer pair No 6 (Table SI) which was designed 
based on the 5' flanking sequence and the conserved 3' flanking 
sequence specific to the candidate SalFAD3.LA2 gene (Figure S3). 
Comparison of the coding and genomic DNA sequences indicated 
that the candidate SalFAD3.LA2 gene also contained 8 exons and 7 
introns (Figure 5). Variation in the length of the introns was 
observed between LA 2 and la 2 (Figure 5). For instance, the third 
intron of LA was 530 bp, while that of la was 1165 bp. 

Sequence alignment of LA 1 and LA 2 indicated that LA 1 
harboured a 9 bp deletion at position 46 (Figure 5 and Figure 
SI), which resulted in the loss of the three amino acids glycine- 
arginine-lysine at position 16. In addition, 77 point mutations were 
observed between LA and LA (Figure SI), of which 19 mutations 
led to amino acid changes (Figure S5). The candidate SalFA- 
D3.LA1 and SalFAD3.LA2 genes exhibited differences in the cloned 
5' flanking sequences (Figure S2), but had the same 3' flanking 
sequences. Variation in the length of the introns was observed 
among the four alleles LA 1 , la 1 , LA and la (Figure 5). 



Phylogenetic analysis based on the polypeptide sequences 
encoded by LA 1 and LA 2 implied that SalFAD3.LAl and SalFAD3. 
LA2 genes in yellow mustard were clustered with FADS genes in 
Brassica species (Figure 6). The SalFAD3.LAl gene was grouped 
together with the FAD3 genes of B. oleracea (Genbank accession 
No.AGH20189), the C genome in B. napus (BnaC.FADSb, Genbank 
accession No.AFJ19037.1) and B. juncea (Genbank accession 
No.ADJ58020.1), whereas the SalFAD3. LA2 gene was in the 
same cluster with the FAD3 genes of B. rapa (BraA.FAD3a, BRAD 
accession No. BraO 18348) and the A genome in B. napus 
(BnaA.FAD3a, Genbank accession No.AFJ19039. 1). 

Co-segregation of the SalFAD3.LA1 and SalFAD3.LA2 
Allele-specific Markers with C18:3 Contents in the F 2 
Populations 

Primer pair 12 (Table SI) produced co-dominant markers of 
742 bp, 510 bp, 626 bp and 1273 bp specific for LA 1 , la 1 , LA 2 and 
la 2 , respectively, which co-segregated with the C18:3 content in all 
of the F 2 populations. In the cross of Y1127 (low)xY1035 (high), 
all of the nine possible genotypes were identified using the markers 
specific for LA , la , LA and la (Figure 7 A; Table 3). The 
homozygous F 2 plants (LA 1 LA 1 LA 2 LA 2 ) had a significantly higher 
CI 8:3 content (average: 17.1%) than the heterozygous F 2 plants of 
LA 1 la 1 LA 2 LA 2 (average: 13.1%) (t= 6.12, p<QM) and of LA 1 LA 1 - 
LA 2 la 2 (average: 15.5%) (t = 2.23, /> = 0.04) (Table 3). The 
homozygous F 2 plants of LA 1 LA 1 la la had a higher average 
CI 8:3 content (14.0%) than those of la 1 la 1 LA 2 LA 2 (average: 9.1%) 
(t = 5.81, /K0.01). In the cross of Y1127 (low)xY514 (medium), 
the three genotypes for C18:3 content were differentiated with the 
markers spec ific for LA 1 and la (Figure 7B; Table 3). The average 
CI 8:3 content of the homozygous F 2 plants (LA LA la la) was 
12.7%, which was significantly higher than the heterozygous F 2 
plants (LA 1 la 1 la 2 la 2 , average: 9.2%) (t = 5.02, /K0.01) (Table 3). In 
the cross ofY514 (medium) xYl 035 (high), the markers specific for 
LA 2 and la 2 distinguished the three C18:3 genotypes (Figure 7C; 
Table 3). The homozygous F 2 plants (LA 1 LA 1 LA 2 LA 2 ) had an 
average C18:3 content of 15.8%, which was higher than the 
heterozygous F 2 plants (LA 1 LA 1 LA 2 la 2 , average: 13.8%) (t = 2.23, 
/> = 0.04) (Table 3). The SalFAD3.LAl and SalFAD3.LA2 genes co- 
localized with the QTL peaks on Sal02 and Sal 10, respectively 
(Figure 4). A new band was observed in the Fi and F 2 plants with 
the heterozygote genotype of LA 2 la 2 (Figs. 7A and 7C). 
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Figure 7. Identification of different C18:3 genotypes based on the SalFAD3.LA1 and SalFAD3.LA2 allele-specific markers in the F 2 
populations of the three crosses Y1127xY1035 (A), Y1127xY514 (B) and Y514xY1035 (C). M: DNA ladder. Y1127: Low C18:3 line 
(la'la'lcfla 2 ). Y1035: high C18:3 line {LA' LA 1 LA 2 LA 2 ). A. F, (Y1 127xY1035): LA'la'LA 2 ^. Lane 1: la'la'lcfla 2 ; Lane 2: LA' LA' I a 2 la 2 ; Lane 3: LA' LA' LA 2 LA 2 ; 
Lane 4: LA'la'la 2 la 2 ; Lane 5: LA'LA'LA^a 2 ; Lane 6: LA' la' LA 2 la 2 ; Lane 7: la'la'LA 2 LA 2 ; Lane 8: LA' la' LA 2 LA 2 ; Lane 9: la' la' LA 2 la 2 . B. F, (Y1127xY1035): 
LA'la'la 2 W; Lane \: la' la' la 2 la 2 ; Lane 2: LA'LA'la 2 la 2 ; Lane 3: LA'la'la 2 la 2 . C. F, (Y514xY1035): LA' LA 1 LA 2 la 2 ; Lane 1: LA'LA'lc^la 2 ; Lane 2: lVm'L4 2 M 2 ; 
Lane 3: L4'L4'L4 2 fcr\ * indicated new band. The different C18:3 genotypes were identified using primer No 12 (SalFAD3 F and SalFAD3 R) (Table S1) 
designed based on the variation in intron 3 of the SalFAD3 alleles LA', la', LA 2 and la 2 . 
doi:1 0.1 371 /journal.pone.0097430.g007 



Discussion 

The present paper reported on the inheritance and QTL 
mapping of CI 8: 3 content as well as molecular characterization of 
the FAD3 genes in yellow mustard. Linolenic acid content was 
controlled by the nuclear genotype of the embryo in yellow 
mustard as reported in B. napus [8]. Two nuclear gene loci were 
detected and functioned independently and additively to deter- 
mine the total C18:3 content in the seeds. However, maternal 
effects on the CI 8: 3 content couldn't be ruled out since 
appropriate progeny tests were not performed in the present 
study. QTL analysis further revealed that the two gene loci 
SalFAB3.LAl and SalFAD3.LA2 had a different magnitude of effect 
and together explained 96.4% of the total variation for CI 8:3 
content. The residual 3.6% variation of C18:3 content beyond the 
two QTLs could be resulted from maternal and environmental 
effects. It has been reported that temperature, maternal genotype 
and cytoplasm have effects on C18:3 content in B. napus [8-10]. 
The duplication of the FAD3 gene provides additional evidence 
that yellow mustard is a secondary polyploid species as revealed by 
molecular studies [29,30]. The two linkage groups Sal02 and SallO 
containing the SalFAD3.LAl and SalFAD3.LA2 genes, did not share 
any common ILP markers, suggesting the occurrence of extensive 
genomic changes during the speciation of yellow mustard. 

Molecular cloning and sequencing indicated that the SalFA- 
D3.LA1 and SalFAD3.LA2 genes contained 8 exons and 7 introns 
in yellow mustard, which is in agreement with that in B. napus [12] 
and A. thaliana (Locus: AT2G29980, TAIR) [31]. However, the 
molecular mechanism underlying the naturally occurring CI 8:3 
variant in yellow mustard was different from that of the EMS- 
induced C18:3 variant in B. napus and B. oleracea. The FAD3 gene 
with reduced C18:3 content resulted from SNP mutations in B. 
napus [11,12] and B. oleracea [32]. However, the recessive allele la 1 
of the SalFAD3.LAl gene was a loss-of-function mutant due to an 
insertion of 584 bp in exon 3. The inserted fragment contained a 
typical LTR retroelement insertion site (5'-AGAAC-3') at each 
end, suggesting that the inserted fragment might be a remnant of a 
transposable element which had undergone a deletion following 



the insertion event. The recessive aUele laf of the SalFAD3.LA2 
gene was functional and had a CDS encoding the same 
polypeptide sequence when compared with the dominant allele 
LA 2 . However, la 2 was different in intron sequence. It remains to 
be investigated why LA 2 and la 2 controlled a different CI 8:3 
content. The SalFAD3.LAl and SalFAD3.LA2 allele-specific mark- 
ers proved to be useful for identification of different CI 8:3 
genotypes in the present study. 

The phylogenetic analysis based on the polypeptide sequences 
indicated that the LA 1 and LA 2 genes in yellow mustard were 
clustered with the FAD3 genes in Brassica species and A. thaliana. 
Interestingly, LA 1 and LA 2 were clustered into different groups. 
LA was grouped together with the FADS genes of B. oleracea and 
the C genome in B. napus, whereas LA 2 was in the same cluster 
with the FAD3 gene of B. rapa and the A genome in B. napus. In our 
study, the LA 1 gene controlled a higher CI 8: 3 content than the 
LA 2 gene. It was reported that the FAD3 gene of the C genome in 
B. napus also contributed more to the total CI 8: 3 content than that 
of the A genome [11,12]. This suggested that the molecular 
divergence of the LA 1 and LA 2 genes occurred before the 
speciation of yellow mustard and Brassica species. 

In conclusion, our study revealed the existence of two FADS 
gene loci contributing to the genetic variation of linolenic acid 
content in yellow mustard. The SalFAD3.LAl gene was located in 
the linkage group Sal02, while the SalFAD3.LA2 gene in SallO. We 
have cloned the SalFAD3.LAl and SalFAD3.LA2 genes and 
developed allele-specific markers for the detection of desirable 
genotypes, which will be valuable for marker assisted breeding in 
yellow mustard. 

Supporting Information 

Figure SI Alignment of the coding DNA sequences of the 
SalFAD3.LAl alleles LA 1 , la 1 and the SalFAD3.LA2 alleles 
LA 2 and la 2 in yellow mustard. The nucleotide sequence 
alignment was carried out using ClustalW2 (http://www.ebi.ac. 
uk/Tools/ msa/ clustalW2/). 
(PDF) 



PLOS ONE | www.plosone.org 



9 



May 2014 | Volume 9 | Issue 5 | e97430 



FAD3 Genes and Allele-Specific Markers in Yellow Mustard 



Figure S2 Alignment of the 5' upstream sequences of 
the SalFAD3.LAl and SalFAD3.LA2 genes in yellow 
mustard. UP-LA 1 and UP-la represented the 5' upstream 
sequences of the alleles LA 1 and la 1 of the SalFAD3.LAl gene. UP- 
LA and UP-la indicated the 5 ' upstream sequences of the alleles 
LA 2 and la 2 of the SalFAD3.LA2 gene. The nucleotide sequence 
alignment was carried out using ClustalW2 (http://www.ebi.ac. 
uk/Tools/ msa/ clustalw2/). 
(PDF) 

Figure S3 PCR amplification of the SalFAD3.LAl and 
SalFAD3.LA2 genes. A. PCR amplification of the genomic 
DNA sequences of Sa.lFAD3.LAl gene using the primer pair No 5 
(Table SI). Lanes 1-2: 4268 bp fragment of la 1 from Yl 127; Lanes 
3-6: 4534 bp fragment of LA 1 from Y1035. B. PCR amplification 
of the genomic DNA sequences of SalFAD3.LA2 gene using the 
primer pair No 6 (Table SI). Lanes 1-4: 4688 bp fragment of la 2 
from Yl 127; Lanes 5-6: 4042 bp fragment of LA 2 from Y1035. 
(PDF) 

Figure S4 Nucleotide sequences of intron 3 and its 
flanking of SalFAD3.LAl and SalFAD3.LA2 genes in 
yellow mustard. The sequence of exon 3 was underlined in 
red while that of exon 4 was lined in blue. The nucleotide 
sequence of the inserted fragment in exon 3 of la 1 was underlined 
in pink. The new intron splicing site GT in the inserted fragment 
was indicated in green rectangle box. The nucleotides in blue 
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